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(— I I Abstract 

The generation of the single spin asymmetries by the twist 3 effects in 
^ . QCD is discussed. While the short-distance subprocesses is calculable 

H ' in the Born approximation, some properties of large distance partonic 

correlations may be extracted by means of sum rules relating them 
to the DIS spin structure function 52- The simple model for twist 3 
part of the latter is proposed. The numerical analysis seems to favor 
the hypothesis about small value of the gluonic poles contribution 
to the single asymmetries, while the fermionic poles may lead to the 
measurable effects. 



1 Introduction 



The single transverse spin asymmetries are known to be one of the most 
subtle effects in QCD. They should be proportional to mass scale, and the 
only scale in "naive" perturbative QCD is that of the current quark mass. 
The additional suppression comes from the fact, that single asymmetries 
are related to the antisymmetric part of the density matrix. Due to its 
hermiticity, the imaginary part of scattering amplitude is relevant. As a 
result, the spin-dependent contribution to the hard scattering cross section 
starts at the one-loop level only. More exactly, it is due to the interference 
of the one-loop spin-flip amplitude and Born non-flip one. At the same time, 
the Born graphs provide a leading approximation to the spin-averaged cross 
section and the asymmetry is proportional to 

The more accurate application of perturbative QCD, including twist-3 
effects, results in a completely different picture. The twist-3 quark-gluon 
correlations give rise to the QCD single transverse asymmetries suppressed 
neither by the quark mass (it is substituted by the hadronic one M) nor by 
as (see e.g. 0, |^ and refs. therein). The collective gluon field of the polarized 
hadron, in which quark is propagating, provides the latter by the mass of 
order that of the hadron. This field is also the source of the phase shift and 
the loop integration in the short- distance subprocess is no more required. 

2 QCD Factorization for single asymmetries 

This qualitative picture is based 0] on a self-consistent approach to the single 
asymmetries in twist-3 QCD. As a result, the parton-like expression was ob- 
tained. A short- distance part is calculable in perturbative QCD with slightly 
modified Feynman rules. The imaginary part is produced in the Born ap- 
proximation due to the extra light-cone integration emerging at the twist-3 
level. A long-distance contribution is described by new two-argument parton 
matrix elements, the so-called quark-gluon correlators (or correlations). The 
latter should, in principle, be determined experimentally from a "partonome- 
ter" process, just like the ordinary parton distributions are determined from 
the deep inelastic scattering. It is important that a mass parameter implied 
by the correlators is of order of the polarized hadron mass. 

Consider a hard inclusive process with a transverse polarized nucleon. 
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The term in the cross section proportional to twist-3 correlators can be ex- 
pressed in the form 

das = J dxidx2^Sp[S^{xi,X2)Tf_,{xi,X2)], (1) 

where Sf^{xi,X2) is the coefficient function of parton subprocesses with two 
quark and one gluon legs; T^(xi,X2) depends on parton correlators: 

T^ixi, X2) = ^{pi1^s^bA{xi,X2) - i7pe™ X2)), (2) 

where e''^*^^ = ^''^"'^ SaPip-, is the covariant hadron polarization vector and 
M is the hadron mass. Two-argument distributions (correlators) 



M J ZTT 

are real and dimensionless. They possess symmetry properties which follow 
from T-invariance 

bA{Xl,X2) = bA{x2,Xi), bv{Xi,X2) = -bv{x2,xi). (5) 

It is convenient to decompose these functions to the singular and regular 
pieces: 



bA{xi, X2) = ipA{xi)6{xi - X2) + b\{x2, xi), (6) 

Xi - X2 

The different form of the singular pieces is dictated by different symmetry 
properties and is of different physical origin. The "short-range" (5-correlation 
in bA is due to the ordinary derivative coming from the transverse momentum. 

At the same time, the pole residue in by 

(pv{xi) = {Xi - X2)bv{xi,X2)\xi=X2 (7) 
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by use of the equation of motion is just the matrix element appearing in the 
"gluonic pole" contribution of J. Qiu and G. Sterman [Q. These authors 
suggested such a term because of the pole of the gluonic propagator in the 
partonic subprocess, producing the imaginary phase which is necessary to 
have the single asymmetry. However, it is natural to use this term for the 
pole of the correlation as well. 



3 Sum rules for partonic correlations 

The quark equation of motion for each flavour allows to relate them to the 
"transverse" spin-dependent quark distribution 

^ J dxdy{bA{y,x)[a{x) + cT{y)] - [cr{x) - a {y)]bv {y , x)) = (8) 

2 j dxa{x)xc^{x). 

Here a{x) is an arbitrary test function. The independence on the choice 
of the vector n, fixing gauge and transverse direction results in the relation: 

^ / dxdy{bA{y,x)^^^^^^ = I dxa{x){ci{x) - 4{x)), (9) 
27r J X — y J 

where cf (x) is the most familiar distribution of the longitudinally polarized 
quarks. All the integrals in (PJ^) are performed in the regions \x,y,x — y\ < 
1. For ordinary parton distributions positive arguments correspond to the 
quarks and the negative ones - to the antiquarks. In the Born approximation 
the structure functions gi and g2 are: 

gi{x)+g2{x) = ^^e2(cT(x) +ct(-x)) (10) 
^ / 

^i(^) = 9 II e^icLix) + cl(-x)) 
^ / 

The quark-gluon correlations in our approach combine the contributions, 
related to the terms of twist 2 and 3 in the operator product expansion. 
However, it appears possible to separate them just from the physical reasons. 
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Note that matrix elements, containing the covariant derivative, are actually 
not gauge invariant. This is because the derivative and the quark field it is 
acting on are taken in the different points at the light cone. One can easily 
pass to the gauge-invariant form by shifting gluon field to the point of quark 
field and expressing the compensating term A^{Xi) — A^{X2) in terms of the 
gluon strength (the latter is posiible because the axial gauge is used). This 
contradicts the earlier statements that transverse momentum and gluon 
field are combined together in a gauge-invariant way. One may separate 
these contributions, extracting the short- distance part of the correlation (0), 
coming from the transverse momentum and the part of gluon field composing 
the covariant derivative. In this approximation 

bA{xi,X2) = ipA{xi)6{xi - X2), bv{xi,X2) = 0. (11) 

The spin structure functions are then completely determined by the trans- 
verse momentum distribution 



9i{x) + g2{x) = - E ^ (12) 

Z ^ X 

which may be excluded from these equations. As a result, one get the simple 
differential equation, whose obvious solution 

9iix)+g2{x)= f'dx^ (13) 

Jx X 



after calculation of the moments 

/ g,(x)+g2{x))dx = 0- (14) 

JO n + 1 

is nothing else than Wandzura-Wilzcek (WW) |^ sum rules. Originally they 
were derived in the framework of the operator product expansion for the 
"twist-2" part of (72 related to the totally symmetric operators. Therefore, 
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the twist-2 approximation corresponds to the taking into account the trans- 
verse degrees of freedom of polarized quark (the same resuh was obtained 
recently by P. Mulders and R. Tangerman using the different approach |^). 

One may easily calculate the gluonic corrections to these sum rule [|], 
defining the twist-3 part ^2, combining (H), (||), and (0) and taking into 
account, that only regular piece of 6^ contributes: 



1 n 

x''g2{x)dx= x"( — —gi{x) + g2{x))dx = (15) 
JO n + I 

— — — / dxidx2^e)[-hv{xi,X2){x1~^ - x^'^) + 

7r[n + 1) J\xi,x2,xi-X2\<i s 2 



f 



bAixi,X2)(l)nixi,X2)], 0„(a;,?/) = -(x" -y" ), n = 0,2... 

X — y I 

Let us consider the important particular cases. For the analysis of the 
famous Burkhardt-Cottingham |0 sum rule in = 0) one may use the ex- 



pression for g2, generated by the (regularized) 5-function as a test function 



m 



f \ 1 /■ , , b''{x,y)+b''{y,x) 
g2{x) = — / dxdy . (16) 

2tT J\x,y,x--y\<l X — y 

One of the possibilities for BC sum rule violation is the long-range singu- 
larity 6{x) 0.1f g2 contains such a term proportional to S{x), which can never 
be observed experimentally, it will give a non-zero contribution to the inte- 
gral and therefore will violate the BC sum rule. Note that such a situation 
was first noticed by Ahmed and Ross in their pioneer paper on spin effects in 
QCD |T^, where they also mentioned a similarity with the Schwinger term 
sum rule for the longitudinal structure function. However, to obtain such a 
behavior, one also needs to have a singularity in e.g. 

b^ix, y) = 6{x)(t){x, y) + 6{y)(j){y, x), (17) 

where is regular function. Such a singularity should result in meaningless 
infinite fermionic poles (when x = 0) contribution to single asymmetries. 

While these arguments allow one to impose some restrictions for g2 start- 
ing from single asymmetries, the most interesting for our purposes is the 
'opposite' problem: what can we learn about twist 3 single asymmetries 
having some information about twist 3 contribution to g2- 
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4 Quantitative estimates of single asymme- 
tries and model for twist-3 part of §2 



The first numerical estimate was performed a couple of years ago [12]. The 



starting point was the value of dxx"^ g2{x) , calculated in the framework of 



QCD sum rule method |T3[. Using the model assumptions about the con- 
stance of chromomagnetic field inside the nucleon, the authors estimated 
quantitatively the gluonic poles contribution to the single asymmetry of di- 
rect photon production. 



To understand this problem better it instructive to consider (JTBI) for n=2 



1 

x'^g2{x)dx = (18) 
/ dxidx2^e'jbv{xi,X2){xi - X2). 

OTT J \xi,X2,xi—X2\<l J 

Note that only vector correlation contribute to this first nonzero moment. 
This fact, known for many years, is especially clear from the non-local OPE 



14] in the co-ordinate space: 



{T7/>7'^V'(x)V'7X2/)}[H = e^'""'{id;D\x, y)i,{x)^^^,U{x, y)^(y) (19) 

- gD%x,y)i>{x)-f^j5 f dt{t - l){x - y)^U{x, z)Fp^{z)U{z,y)^{y) 

Jo I 

- igD'{x,y)^P{x)j^J^ dt^U{x, z)Fpx{z{t)){x - y)^U{z,y)ij{y)}. 

The term ip{0)'y„'j^Fpx{0)'il!{0) contributing to the second moment of g2 
contains the factor dt{t — |) = 

The sum rule (^) is especially simple when one neglect the regular piece 
of hv'- 

I x^g2{x)dx = -^ I dxidx2^e^f(pv{xi)- (20) 

Performing the integration over X2 one get: 
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/ x^g2{x)dx = -— rfa;i Ve^.v9v/(xi)(2 - (21) 
JO An J-i J 

The contribution of the first term in the bracket coincides, up to a con- 
stant of order 1, with the sum rule derived |]12| by use of the model assumption 
about the constant collective gluon field. This is quite natural, because the 
gluonic pole in the correlation is produced by the integration over large dis- 
tances at the light cone, corresponding to the quasiclassical collective fields. 
The second term, effectively increasing the correlation for given value of g2, 
is produced by the shape of integration region. The latter is related to the 
analytical properties of the correlations and is therefore a pure quantum 
effect. 

The second moment of g2 was calculated, as it was mentioned above, 
in the framework of QCD sum rules|]13[. These calculation was checked 



recently |jT3[ by use of the another interpolating nucleon currents, explicitly 
accounting for the gluonic degrees of freedom. The results are compatible: 
this moment is small (about 10~^) and negative for the proton case, while for 
the neutron it is about of order of magnitude larger. Consequently, the single 
asymmetry of direct photon production is claimed to be small and strongly 
flavour- dependent . 

One should note, however, that only gluonic poles contribution calculated 
by J. Qiu and G. Sterman was taken into account. At the same time, the 
contribution of fermionic poles (when imaginary part is produced by the pole 
of quark propagator) is related to the correlation La, which is not manifested 
at all for the second moment of (72 • 

While the fermionic poles were suggested originally as a source of single 
asymmetries 0, the gluonic poles were discovered later ||^ and claimed the 
dominant contribution (because it is more easy to emit the zero-momentum 
gluon than zero momentum quark; note, however, that such emission is 
known to produce the IR divergencies and may therefore cancel somehow). 

Now, when QCD sum rules tell us that gluonic poles are, contrary to 
these expectations, small, there are two opportunities. First (seems to be 
now silently accepted) is to believe that the contributions of fermionic poles 
are even smaller. Second, which we study here, is to come back to fermionic 
poles and investigate them more seriously. 
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To do this one may, say, repeat the QCD SR calculation substituting F — >■ 
F, 7o- — i> 7cr75. This operator, as it was already explained, does not contribute 
to g2- However, its non-local generalization is directly related to h^. If one 
observe its magnitude larger enough than that of vector operator, this would 
be the strong support that only the gluonic poles are small. The same may 



be probably performed by modifying the existing lattice ||T6[ calculation of 
the second moment of g2- 

As we do not know the value of this matrix element at the moment, let 
us try to perform at least the model estimates. To do this, another sum rule, 
relating the first moments of gi and (72, is extremely useful. It is known for 
several years 0, although the role of valence quarks was emphasized only 
recently l]!^ . 



The applicability of the correlations to the single asymmetries is possible 
because they are the (generalized) distribution functions on the light cone, 
accumulating the large-distance dynamics. This allows one to apply (^J^) 
with the odd test function and produce sum rules, which can never appear in 
the framework of the operator product expansion, because they corresponds 
to the "wrong" moments (the positive parity of two photons pick up the 
even combinations only). In these cases the combinations cl,t{x) — cl,t{—x) 
appear. They are just the valence quark contributions. As a result, the sum 
rules (|T5|) are valid for n = 1,3... with the gi,g2 in the l.h.s. being replaced 
by their valence parts. The case of n = 1 is especially interesting, because 
the contribution of the correlations in the r.h.s. appears to be equal to zero. 

dxx{g'(''\x) + 2^™^(x)) = (22) 

Jo 

This sum rule allows one to relate such a "model" entities as valence 
quark starting from the QCD correlations. The situation here is in some 
sense opposite to, say, Gottfried sum rule. Instead of "model" derivation for 
measurable functions, we have (more or less) rigorous QCD derivation for 
"model" functions. Note that this sum rule is not affected by the possible 
singularity, spoiling the validity of BC sum rule . 

The valence contribution to g2 may be measured, in principle, by consid- 
ering the semi-inclusive asymmetries in DIS on transversely polarized target, 
in complete analogy with the studies of longitudinal valence distributions at 
CERN |T^. Although the transverse polarization is, generally speaking, the 
non-probabilistic twist 3 effect, the parton-like formula ( |T0| ) for DIS allows 
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one to extract the valence parts of ct (and, by use of (H), of the correlations), 
applying the same formulas as for the longitudinal polarization. One should 
just change the longitudinal semi- inclusive asymmetries for the transverse 
ones, and cl to ct- 

The alternative derivation of this sum rule 0] leads to its validity for total 
structure functions. One may guess that it is valid for quark and antiquarks 
separately. The other opportunity is [20| that the quark-gluon correlations 
are negligible, if the signs of xi and X2 are opposite (i.e., the correlation of 
gluon and quark-antiquark pair is considered). It is interesting, that in this 
case one should change 2 — |xi| — 1 in the sum rule for the second moment 



0]) and get the expression completely analogous to its model derivation [|12 



Anyway, let us neglect in the following the difference between valence and 
total contribution to (?2- 

As a result, the WW sum rules appear to be exact (in the leading twist 3 
approximation) for the two moments. Consequently, the deviation of g2 from 
its WW value (twist 3 piece) has these moments equal to zero. This function 
should change sign at least two times inside the region (0,1). The simi- 
lar behavior was predicted for the singlet piece of the tensor spin structure 



function 21 



Let us incorporate this property to the following parametrization of the 
twist 3 contribution, where, as usual, the Regge behavior for small x and the 
phase space for x ~ 1 are manifested: 



g2{x) = Cx'^ix - ai) (a; - 02) (1 - xf (23) 
Imposing the conditions 

/ dxg2{x) = / dxxg2(x) = (24) 
JO Jo 

allows us to find the crossing points ai and 02 in terms of a and jS. At the 
final stage, we shall calculate the second moment of g2- 

We performed these calculation for (3 = 4 and a changing in the interval 
[0,1). As one should not take this parametrization too seriously, the general 
features, not very sensitive to the actual values of parameters, are most 
important: 

i) The function has the negative minimum at a; ~ 0.12 ± 0.05 and the 
positive maximum at a; ~ 0.6 ± 0.05. 
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ii) The second moment is rather small and negative - about — 10~^C. 

iii) The ratio of the value of function in the minimum to the value of 
the second moment is changing in the interval —(70 — 130), while for the 
maximum it is +(10 — 20). 

One may expect, that the small value of the second moment, predicted 
by QCD sum rules, is reflecting the strong oscillations rather than small nor- 
malization. This is also compatible with the dominance of axial correlation, 
absent for the second moment but manifested in another moments and in 
the fermionic poles contributions to single asymmetries. 

Taking into account the negative QCD SR value, one should expect the 
maximal value of ^2 at x ~ 0.1 to be positive. 



^^(x = 0.12 ±0.05) = 0.16 ±0.10; (25) 
g^{x = 0.12 ±0.05) = 0.8 ±0.4; 

While the preliminary data of E143 experiment reported in Morionde were 
negative in this region (supporting the positive value of the second moment), 
the more recent data [E^ are indicating (for one of two spectrometers) the 



positive values compatible with (25). However, they are still compatible with 



zero as well, and the more precise HERMES data are most welcomed. 



5 Short-distance subprocesses 

Let us discuss briefly the calculation of the short-distance QCD subprocess. 
It is determined by Born graph (the imaginary part is produced by the pole 
of the propagator while the integration over light-cone variable is performed) 
and the computation is straightforward. Note that a number of fermionic 
poles contributions to quark-gluon subprocesses is already calculated . 
This offer a possibility not to limit the future experimental studies by the 
asymmetries of the direct photon production. 

I would like to stress the special role of the single asymmetries of jet 
production. This asymmetry is just 

A = ^^1^, (26) 
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where ctl (o"_r) is just the cross-section for jets produced to the left (right) 
of the normal to scattering plane. This does not require the full analysis of 
jet structure, like for handedness studies [^], and only the trust direction is 
necessary. At the same time, the cross-sections are much larger than in the 
photon case. In the proposed experiment HERA — N, when 800 GeV proton 
beam would scatter on the polarized nucleon target, the expected statistical 
error for the jet left-right asymmetry is less than 10~^. This would allow one 
to measure the asymmetry of order 10~^, expected for "pessimistic" estimates 
based on the values (^), 

From the theoretical point of view, these effects are just the asymmetries 
of quark and gluon production. They are insensitive to the fragmentation 
effects and provide as clear test of QCD, as a photon asymmetry. 

Another interesting process, for which calculations are already performed, 
is the dilepton left-right asymmetry. The dilepton mass variation may be 
used to probe the quark-gluon correlation in various regions of the ligh-cone 
momentum fractions. At the same time, the background problem can be 
solved more easy in this case. 

The quark gluon subprocesses contribute also to the single asymmetries 
of pion production. It is interesting, that the well-known problem of flavour 
dependence can be naturally solved, if the contribution of the fermionic poles 
is much larger than that of the gluonic poles. 

The non-existence of the gluonic poles leads to the significant numerical 
difference between the single asymmetries of the gluon 



_ bA{0,x)-bv{0,x)MpT 
fix) m\ 

{Cf - CAim - ^f) ((4 + 1)Ca/2 - Cpxp{\ - xpf) 
^ Cp{\\x\) {xpiCp-CAm-iX^xDCpliy 

and quark 
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A, 



bA{0,x) -bv{0,x) MpT 
" fix) m2 



T 

1 

X- 



Cf{Cf{xUI - xf) + 4/2) + Ca/2(2 - 4xf + 3x|. - 4)) 
- 1)) - C^C^/2(4 - - + 

+Cy4(x|, - 5xp + - 14x| + 9xf - 4) 



X 



production on the transverse polarized nucleon by the gluon from the unpo- 
larized one. Namely, the latter is several times larger [p5|. One may expect, 
that the pion single asymmetry is due to the pions resulting from the quark 
(not gluon!) fragmentation. As a result, the tt"*" meson asymmetry is related 
to the correlation of gluon and u-quark, while tt~ meson asymmetry is related 
to the correlation of the gluon and rf-quark; 7r° can be produced by u- and 
(i-quarks with equal probability. 

Making use of sum rules relating the correlations to the "ordinary" quark 
distributions (H), valid for each flavour separately, one can easily get 

Am Ad Au + Ad 
~ , A^- ~ — — , A^o ~ r~^- 

u d u + d 

Note that mirror asymmetries for At^+ and A^- and A^^o = 1/3^4^+ can 
be obtained if u = 2d and Au = —2Ad, what is not far from experimental 
observation. 



6 Conclusions 

The earlier quantitative estimates of the twist-3 single asymmetries are deal- 
ing, in fact, only with the contributions of the gluonic poles, their smallness 
therefore should be interpreted as the argument against the original hypoth- 
esis of J. Qiu and G. Sterman about the domination of gluonic poles. The 
relative strength of fermionic and gluonic poles may be checked by the calcu- 
lation of "dual" quark-gluon matrix element using QCD sum rules or lattice 
simulations. 

As the solution of this problem by use of nonperturbative methods was 
not found yet, the models of twist 3 effects are especially important. We 
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are suggesting the simple model for the twist 3 part of g2 structure function. 
It incorporates the Regge and phase space behavior, Burkhardt-Cottingham 
sum rule as well as another sum rule relating the first moments of gi and 
g2. It is found, that the smallness of the second moment is the result of 
the strong oscillations required by sum rules, while its absolute value may 
be fairly large at the region of ,t ~ 0.1. The values of twist-3 part of g2 in 
this region are predicted in the proton (about 0.15) and neutron (about 0.8) 
cases, using the QCD sum rules calculations for the second moment. Such a 
behavior is supporting the domination of the fermionic poles over the gluonic 
ones, because it is the second moment, where the correlation producing the 
fermionic poles does not contribute. 

The existing perturbative calculations allow one to study the various sin- 
gle asymmetries in the proton-nucleon scattering at 800 GeV (proposed ex- 

— * 

periment HERA — N) . The left- right asymmetries of jets and dileptons are 
especially interesting, as providing good statistics and offering the opportu- 
nity of clear tests of twist 3 QCD. 
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